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ABSTRACT 



X 

O Context. The ~ 1700 year old PSR B0540— 69 in the Large Magellanic Clouds (LMC) is considered the twin of the Crab pulsar 
■ because of its similar spin parameters, magnetic field, and energetics. PSR B0540— 69 (V ~ 22.5) is also one of the very few 
C/5 pulsars for which both optical pulsations and polarised emission have been measured. Its optical spectrum is fit by a power-law, 
, ^ , , ascribed to synchrotron radiation, like for the young Crab and Vela pulsars. At variance with them, however, a double break is 

required to join the X-ray and optical power-law spectra, with the first one possibly occurring in the near ultraviolet (nUV). 
t— h Aims. Near-infrared (nIR) observations, never performed for PSR B0540— 69, are crucial to determine whether the optical 
. power-law spectrum extends to longer wavelengths or a new break occurs, like it happens for both the Crab and Vela pulsars 
' in the mIR, hinting at an even more complex particle energy and density distribution in the pulsar magnetosphere. 
I/"} ' Methods. We observed PSR B0540— 69 in the J, H, and Ks bands with the Very Large Telescope (VLT) to detect it, for the first 
time, in the nIR and characterise its optical-to-nIR spectrum. To disentangle the pulsar emission from that of its pulsar wind 
' nebula (PWN), we obtained high-spatial resolution adaptive optics images with the N AOS-CONIC A instrument (NACO). 
, Results. We could clearly identify PSR B0540— 69 in our J, H, and Ks-band images and measure its flux (J=20.14, H=19.33, 
1 Ks=18.55, with an overall error of ± 0.1 magnitudes in each band). The joint fit to the available optical and nIR photometry 
with a power-law spectrum F v tx v a gives a spectral index a = 0.70 ± 0.04, slightly more precise than measured in the optical 
only. This suggests that there is no spectral break between the optical and the nIR. We also detected, for the first time, the PSR 
B0540— 69 PWN in the nIR. The comparison between our NACO images and Hubble Space Telescope (HST) optical ones does 
not reveal any apparent difference in the PWN morphology as a function of wavelength. The PWN optical-to-nIR spectrum is 
p\ . also fit by a single power-law, with spectral index a — 0.56 ± 0.03, slightly flatter than the pulsar's. 
^ ' Conclusions. Using NACO at the VLT, we obtained the first detection of PSR B0540-69 and its PWN in the nIR. Due to the 
small angular scale of the PWN (~ 4") only the spatial resolution of the James Webb Space Telescope ( JWST) will make it 
possible to extend the study of the pulsar and PWN spectrum towards the mid-IR. 

Key words. Optical: stars; neutron stars: individual: PSR B0540-69 



1. Introduction 

Historically discovered in radio, pulsars emit across the 
whole electromagnetic spectrum, from the radio, to the 
optical (Mignani 2011a), to the X-rays (Becker 2009), and 
to high-energy 7-rays (Abdo et al. 2010a). The pulsar op- 
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tical emission is due to the combination of two different 
processes: synchrotron radiation from relativistic electrons 
in the neutron star magnetosphere and thermal radiation 
from the cooling neutron star surface. Their signature is 
recognised by a power-law (PL; F v oc v~~ a ) spectrum with 
a i; and by a Rayleigh- Jeans (R-J) continuum, respec- 
tively. The relative weight of these two spectral compo- 
nents depends on the pulsar age. Young (^ 10 kyr old) 
pulsars feature dominating magnetospheric emission, ex- 
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tending to both the near Infrared (nIR) and the near ul- 
traviolet (nUV), while the middle-aged ones (i; 100 kyr 
old) feature both magnetospheric and thermal emission, 
with the former dominating in the nIR and the latter in 
the nUV. 

Very little is know about the pulsar IR emission. Only 
5 out of the 12 pulsars identified in the optical, have been 
also detected in the nIR. These are the young pulsars 
Crab (e.g., Sandberg & Sollerman 2009; Tziamtzis et al. 
2009), Vela (Shibanov et al. 2003), and PSR B1509-58 
(Kaplan & Moon 2006), and the middle-aged pulsars 
PSR B0656+14 and Geminga (Koptsevich et al. 2001). 
Moreover, the Crab and Vela pulsars have been also 
clearly detected in the mIR by Spitzer (Temim et al. 
2009; Danilenko et al. 2011), with a possible marginal in- 
dication of detection also for Geminga (Danilenko et al. 
2011). Spitzer also observed PSR J1124-5916 (Zyuzin et 
al. 2009), PSR J0205+6449 (Slane et al. 2008), and PSR 
J1833-1034 (Zajczyk et al. 2011) but it could only de- 
tect the surrounding pulsar- wind nebulae (PWNe). Spitzer 
also detected the Crab PWN, which is the only one de- 
tected both in the nIR (Sandberg & Sollerman 2009) and 
in the mIR (Temim et al. 2009), together with the PSR 
J1833-1034 PWN (Zajczyk et al. 2011). A nIR source 
has been recently associated with the 38.5 ms X-ray pul- 
sar IGR J18490-0000 (Gotthelf et al. 2011; Ratti et al. 
2010; Curran et al. 2011). However, the identification with 
the pulsar has never been confirmed so far. Detecting more 
pulsars in the IR is crucial to search for continuity or spec- 
tral breaks at wavelengths longer than the optical band 
and, in turn, to constrain the properties of particles in 
the pulsar magnetosphere on a broader spectral range. 
Obviously, young and energetic pulsars are the best tar- 
gets. 

Interestingly enough, the X-ray /radio pulsar PSR 
B0540-69 (Seward et al. 1984; Manchester et al. 1993) 
in the Large Magellanic Clouds (LMC) has not been ob- 
served yet in the IR. With a spin-down age of ~ 1700 yr, 
it is one of the youngest pulsars so far and is similar to 
the Crab in its spin period (P=50 ms), dipolar magnetic 
field (B ~ 5 x 10 12 G), and rotational energy loss rate 
(E ~ 1.5 x 10 38 erg s _1 ). Moreover, it is the second bright- 
est pulsar (V ~ 22.5) identified in the optical (Caraveo et 
al. 1992) and one of the very few for which both pulsa- 
tions (Gradari et al. 2011) and polarisation (Mignani et 
al. 2010) have been measured in this band. High-spatial 
resolution images taken with the Hubble Space Telescope 
(HST) clearly resolved PSR B0540— 69 from its compact 
(~ 4" diameter) PWN (Mignani et al. 2010), allowing the 
most accurate measurement of the pulsar optical spec- 
trum. The optical PL slope (a — 0.70 ± 0.07) is similar to 
the X-ray one (a x = 0.92 ± 0.11; Kaaret et al. 2001), but 
the optical fluxes fall orders of magnitudes below the X-ray 
PL extrapolation, as observed in most pulsars (Mignani 
et al. 2010). Interestingly enough, while for both the Crab 
and Vela pulsars a single break is sufficient to join the op- 
tical and X-ray spectra, two breaks are required for PSR 
B0540— 69. This unique behaviour hints at a complex par- 



ticle energy and density distribution in the pulsar magne- 
tosphere. PSR B0540— 69 has never been observed so far 
in the nIR, while in the mIR the coarse spatial resolu- 
tion achievable with Spitzer made it impossible to resolve 
the pulsar from its bright PWN (Williams et al. 2008). 
Thus, we do not know whether a single PL smoothly joins 
the nIR and optical spectra or the PL slope changes once 
again, like it happens for the Crab and Vela pulsars in 
the mIR (Tziamtzis et al. 2009; Danilenko et al. 2011). 
In the latter case, an even more complex scenario would 
emerge for PSR B0540— 69, with three breaks across its 
nIR-to-X-ray PL spectrum, a challenging one for pulsar 
magnetosphere models. 

Here, we present the results of the first nIR observa- 
tions of PSR B0540— 69, recently performed by our group 
using adaptive optics at the Very Large Telescope ( VLT) . 
This paper is organised as follows: observations, data re- 
duction and analysis are described in Sect. 2, while results 
are presented and discussed in Sect. 3 and 4, respectively. 
Conclusions follow. 

2. Observations and data reduction 

2.1. Observation description 

PSR B0540— 69 was observed in service mode on October 
29th and December 9th 2010 at the ESO Paranal 
Observatory with NAOS-CONICA (NACO), the adap- 
tive optics (AO) imager and spectrometer mounted at 
the VLT Yepun telescope. In order to provide the best 
combination between angular resolution and sensitivity, 
we used the S27 camera with a pixel scale of 0'.'027 
(28" x 28" field-of-view). As a reference for the AO cor- 
rection we used the GSC-2 star S0101001160715 (Lasker 
et al. 2008), located ~ 22"away from our target, with 
the visual (4500 — 10000 A) dichroic element and wave- 
front sensor. We observed PSR B0540— 69 in the stan- 
dard J (A = 1.265,um; AA = 0.25^m), H (A = 1.66^m; 
AA = 0.33^m) and K s (A = 2..8/im; AA = 0.35^m) fil- 
ters. To allow for subtraction of the variable IR sky back- 
ground and avoid possible ghosts from the bright reference 
star (V = 13.05) at the edge of the NACO field-of-view, 
we split each integration in sequences of short randomly 
dithered exposures with detector integration times (DIT) 
of 20 s and 24 jittered exposures (NDIT) along each node 
of the 5-point dithering pattern. This corresponds to a net 
integration time of 2400 s per filter. 

For all observations, the seeing conditions measured 
by the Paranal differential image motion monitor (DIMM) 
were on average below ~ 0"8 and the airmass was ~ 1.4, 
i.e. about the minimum possible due to the low declination 
of our target. Sky conditions were photometric in both 
nights according to the ESO ambient condition monitofl 
Atmospheric water vapour content and wind speed were 
well within tolerance values. Night (twilight flat-fields) 
and daytime calibration frames (darks, lamp flat-fields) 

http: / / archive.eso.org/asm/ambient-server 
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Fig. 1. 15" x 10" J-band image (2400 s) of the region 
around PSR B0540-69 and its PWN taken with NACO 
(top panel). The pulsar is the bright point source at the 
centre of the PWN. The H and Ks-band images are also 
shown in the mid and bottom panels, respectively. 



were taken daily as part of the NACO calibration plan. 
For both nights, standard stars from the Persson et al. 
(1998) fields were observed. 



2.2. Data reduction 

Calibration data were processed using the last version of 
the ESO NACO pipelin|l For each of the three filters 
separately, we dark-subtracted, flat-fielded and cleaned 
for hot pixels the NACO science images with the pro- 
duced master dark and flat field frames. We also applied 
the correction for the odd-even column pattern affecting 
the NACO images. When available, we also used twilight 
flats of similar exposure level as the sky background to 
minimise non-linearity effects and the fixed pattern noise. 
Otherwise, we used a lamp flat. We constructed the sky 
frame via a weighted median stacking of the intermediate 
science products. 

Since the jitter offsets along each node of the 5-point 
dithering pattern are comparable to the extension of the 
SNR around PSR B0540— 69, a contamination was intro- 
duced in the computed sky frame. We corrected this con- 
tamination using the following procedure. We decomposed 
the computed sky frame Stot into several components: the 
constant median value M, a large-scale two-dimensional 
gradient G, a component containing noise N, and a resid- 
ual frame F containing the clean sky S c and the contam- 
ination C produced by the bright SNR, with F = S c + C: 

S tot = M + G + N + F. (1) 

We subtracted the median M and the gradient Gfrom 
the total sky frame Stot and we smoothed the resulting 
frame using cells of 10 pixel size (0"27) to separate the 
residual pixel-to-pixel component N) from the large-scale 
background structure to derive F. For the separation of 
the clean sky S c from the contamination C we replaced 
flux values above a certain threshold by a common fraction 
of their flux: 

( RF : F>=T 
[ F : F < T 

where the threshold T and the flux scaling ratio R were 
free parameters. The best values for T and R have been 
found itcratively to minimise residuals of the contamina- 
tion C and avoid over-correction. Finally, we replaced F 
by S c in Eq. [TJ to reconstruct the corrected sky frame. 

We finally subtracted the cleaned sky as part of the 
co-adding of aligned images and bad pixel correction pro- 
cedure using the NACO pipeline science recipe. The final 
image quality, as measured from the median full-width 
half maximum (FWHM) of all star-like objects in the fi- 
nal product frame, is ~ 0"16, ~ 0"12, and ~ 0"11 for the 
J, H, and Ks-band images, respectively. 

As a reference for the photometric calibration we used 
the zero-points computed by the NACO pipeline, after cor- 
rection for the effects of the atmospheric extinction using 
the most recents coefficients for the Paranal Observatory 
computed by Lombardi et al. (2011). Since no 2MASS 
stars (Skrutskie et al. 2006) are present in the narrow 
NACO/S27 field-of-view, the astrometric re-calibration 

2 www.eso.org/ observing/ dfo / quality /NACO / pipeline 
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of the NACO frames can only be computed using a grid 
of secondary stars selected from our HST images of the 
pulsar field, tied to the International Celestial Reference 
Frame (ICRF) with a ~ 0'.'12 accuracy (ler), as astro- 
metric reference (Mignani et al. 2010). However, since we 
could unambiguously recognise the pulsar in the NACO 
images from relative astrometry with respect to the HST 
ones, we did not deem necessary to pass through this step. 

Of course, the lack of 2MASS stars in the NACO field- 
of-view also means that we have no way to counter-check 
our photometric calibration against an external catalogue. 
To rule out possible mis-calibrations, we thus checked the 
values of the night zero points with those taken on a 30 
day time window centred on the dates of our observations 
and we found perfect consistency, with an rms ^0.07 mag- 
nitudes in all bands. To rule out possible systematics in 
the flux calibration computed by the NACO pipeline, we 
recomputed the night zero points by matching the instru- 
mental magnitudes of the standard stars with those in the 
Persson et al. (1998) catalogue. For each band, we found 
agreement within the rms of the NACO pipeline zero point 
distribution, after accounting for the atmospheric extinc- 
tion correction. Thus, we assume ~ 0.07 magnitudes as 
the uncertainty on our absolute flux calibration, with the 
uncertainties associated with the atmospheric extinction 
correction being negligible. 

3. Results 

The J-band image of the NACO field centred on the PSR 
B0540— 69 position is shown in Fig.[T](top panel). As seen, 
the pulsar is clearly resolved against its bright PWN and 
it is also resolved in the longer- wavelength H and Ks-band 
images. Thus, we detected PSR B0540-69 in the nIR for 
the first time and is now the sixth pulsar detected at these 
wavelengths. 

3.1. Pulsar photometry 

We computed the pulsar magnitudes in the J, H, and Ks 
bands through customised aperture photometry with the 
IRAF task qphot. We used an aperture of 7 pixel radius 
(~ 0'.'19), which is comparable with the median image 
FWHM and small enough to minimise the flux contamina- 
tion from the surrounding, bright PWN. Since the surface 
brightness distribution of the PWN is not homogeneous, 
and it peaks close to the pulsar, the choice of the areas 
used to compute the sky background is crucial for its cor- 
rect characterisation and for not misestimating the pulsar 
flux. We optimised the sky background computation using 
an annulus located at a distance of 25 pixels (0'.'68) from 
the photometry aperture, to minimise the contamination 
from the wings of the pulsar's PSF, and with a variable 
width of 20-30 pixels as a function of the filter, to account 
for wavelength-dependent variations in the surface bright- 
ness of the PWN. We customised the maximum size of the 
background annulus not to include back/foreground stars 
spatially coincident with the PWN. 
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Fig. 2. Extinction-corrected nIR and optical spectral flux 
densities of PSR B0540-69. The values in the J, H, and 
Ks filters have been computed from our new NACO ob- 
servations while the values in the optical filters (labelled 
with numbers) are derived from the HST observations of 
Mignani et al. (2010). 

To the pulsar fluxes measured in each filter, we then 
applied the aperture correction. For each filter, we com- 
puted it from the growth curve of a few bright field stars, 
starting from the same 7 pixel radius photometry aper- 
ture used for the pulsar. We selected those with no ad- 
jacent star within a 1'.'5 radius and detected not further 
than ~ 5" from the pulsar position to minimise aperture 
correction uncertainties due to the variation of the NACO 
PSF at large off-axis angles, produced by the degrade of 
the AO correction. Since the value of the image quality 
varies with the filter (see Sect. 2.1), as the result of the 
more efficient AO correction at longer wavelengths, it was 
crucial to compute the aperture correction per each fil- 
ter separately. To the measured fluxes, we then applied 
the measured zero points and the atmospheric extinction 
correction. The pulsar magnitudes are J=20.14, H=19.33, 
Ks=18.55, with an error of 0.1 magnitudes which accounts 
for the statistical errors (~ 0.05), as well as for the uncer- 
tainties on our absolute flux calibration (~ 0.07) and the 
aperture correction estimate (~ 0.06). 

3.2. The pulsar spectrum 

We used our J, H, Ks-band magnitudes to characterise, 
for the first time, the pulsar flux in the nIR and com- 
pare its nIR spectrum with the extrapolation of the PL 
which best-fit the optical fluxes measured with the HST 
(Mignani et al. 2010). From the measured pulsar mag- 
nitudes we computed the spectral flux densities in each 
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filter using as a reference the unit conversion factors from 
the CIT photometric system (Elias et al. 1982). We then 
corrected the flux densities for the interstellar extinction 
using an E(B — V) — 0.2 (see Mignani et al. 2010 and ref- 
erences therein) and the interstellar extinction coefficients 
of Fitzpatrick (1999). The results are shown in FigJ^l As 
seen, the IR fluxes nicely follow the extrapolation of the 
optical PL spectrum (ao — 0.7 ± 0.07). A joint fit from 
the HST336W to the NACO K s -band fluxes give a spec- 
tral index ceo.niR = 0.70 ± 0.04. For comparison, a fit to 
the nIR data points alone would yield a spectral index 
OiniR = 0.66 ± 0.2, i.e. fully compatible with that mea- 
sured across the whole spectral range. Thus, we found 
no evidence of spectral breaks between the optical and 
nIR, at variance with what is expected between the opti- 
cal and the X-rays. Our derived spectral index is virtually 
identical to that measured by Mignani et al. (2010) using 
only the HST data. However, thanks to the larger wave- 
length coverage, our spectral fit is slightly more accurate, 
although the larger error bars on the NACO photometry 
do not allow to have a more significant improvement. 

Our data obviously do not rule out the presence of a 
spectral break towards longer wavelengths like that found, 
for instance, in the spectrum of the Crab and Vela pulsars 
(Temim et al. 2009; Danilenko et al. 2011). Observations in 
the mIR would be crucial to search for such break but, un- 
fortunately, the coarse angular resolution of Spitzer does 
not allow one to resolve the pulsar emission against that of 
its PWN. Moreover, the pulsar+PWN Spitzer fluxes in the 
IRAC1-4 (3.6, 4.5, 5.8, and 8.0 fj,m) and MIPS1 (24^m) 
channels (Williams et al. 2008) are about an order of mag- 
nitude above the extrapolation of the optical-to-nIR PL 
spectrum of the pulsar, so that they are obviously non con- 
straining. A spectral break at wavelength as long as 5/im 
in the PSR B0540— 69 spectrum can only be detected with 
the Near Infrared Camera (NIRCam) aboard the James 
Webb Space Telescope (JWST), which has an high enough 
spatial resolution required to resolve the pulsar emission 
from that of its PWN. 

3.3. The nebula morphology 

We clearly detected the PSR B0540-69 PWN in the J, H 
and Kg band images (see Fig. [T]) and, it is now the fifth 
PWN detected in either the nIR or mIR. Similarly, as for 
the pulsar, the PWN is detected for the first time in the 
nIR. We studied both the morphology of the PWN in the 
different nIR bands and its nIR spectrum. 

As seen from Fig. [TJ there are no obvious differences 
in the nebula structure seen in different nIR bands. The 
size of the PWN in the nIR is similar to that seen at 
shorter wavelengths, as seen from the comparison with the 
HST/Wide Field Planetary Camera 2 (WFPC2) images of 
Mignani et al. (2010). Thus, no variation in the nebula size 
is observable when going from the optical to the nIR (Fig. 
[3]). Like in the optical (e.g., Serafimovich et al. 2004), we 
can distinguish two main structures in the PWN. The first 
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Fig. 3. Co-added JHK s -band image of the PSR 
B0540-69 PWN with the contours of the 1995 
J ffST/WFPC2 F555W image (Mignani et al. 2010) over- 
laid. Note that the emission knot at ~ 2'.'2 SW of the pul- 
sar, visible in the HST/WFPC2 contours (epoch 1995.87), 
is not detected in our more recent (epoch 2010.94) NACO 
images. 

one is a bright bar, along a position angle PA^ 45° (mea- 
sured North through East), that stretches at both sides of 
the pulsar and is identified as a probable torus (Gotthelf 
& Wang 2000; Serafimovich et al. 2004). Superimposed to 
the bar, south-west (SW) to the pulsar, there is a blob 
of bright emission. This feature is also visible in the op- 
tical HST data (De Luca et al. 2007) and also in the X- 
ray Chandra data (Gotthelf & Wang 2000). However, the 
brightness contrast between the blob and the north-east 
(NE) side of the putative torus seems to be more pro- 
nounced in the X-rays than in the nIR or in the optical. To 
the north of the pulsar, a weak slightly extended emission 
is recognised, also seen in the optical data. Most probably, 
this feature is associated with a possible jet identified in 
the optical by Serafimovich et al. (2004). Like in the op- 
tical, we observe differences in the surface brightness on 
small angular scales (5s l'-'O), which might imply differ- 
ences in the injection of relativistic particles in the PWN, 
possibly due to a different wind geometry. Alternatively, 
they might be produced when Rayleigh- Taylor instabili- 
ties occur at the interface between the PWN bubble and 
the surrounding cold SNR ejecta producing, e.g. filaments 
and/or fingers. We could not detect the bright emission 
knot observed in the optical by the HST ~ 2'.' 2 SW of the 
pulsar (Luca et al. (2007). This can be explained either 
by a sharp break in the spectrum of the knot or, more 
likely, by its intrinsic variability, possibly associated with 
a change in its apparent position. Indeed, both the fading 
of the knot and its apparent displacement are clearly no- 
ticed from the comparison between the multi-epoch HST 
observations of De Luca et al. (2007) and Mignani et al. 
(2010), which span a time base line of ~ 12 years. Possible 
scenarios on both the knot origin and its fading are dis- 
cussed in Lundqvist et al. (2011). 
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Fig. 4. Ks-band image of the PSR B0540-69. The ex- 
traction region used for estimation of the PWN flux in J, 
H, Ks bands and also in the HST bands (Mignani et al. 

2010) is shown with an ellipse. 

3.4. The nebula spectrum 

We used our NACO data to characterise the optical/nIR 
spectrum of the PWN. We measured the JHKs fluxes for 
the whole PWN using an elliptical extraction region of 
size 3'.'4 x 2'.'4 for semi-major and semi-minor axes, respec- 
tively (see Fig. H|) . The contribution from the pulsar and 
stars overlapping the PWN extraction regions was eval- 
uated through aperture photometry and subtracted from 
the total flux. The correction for the airmass and interstel- 
lar extinction has been applied as for the pulsar (see Sectn. 
3.2). Using the same extraction regions as for the JHKs 
bands, we also calculated the corresponding optical fluxes 
using the HST/WFPC2 data of Mignani et al. (2010). As 
before, the contribution from the pulsar and overlapping 
stars were subtracted from the measured values. Following 
the WFPC2 Data Handbook (Bagget et al. 2002) and the 
approach described in Sect. 4.2. of Mignani et al. (2010) 
we corrected the fluxes for the charge transfer efficiency 
losses of the WFPC2 detector. Correction for the inter- 
stellar extinction was applied as usual. The final results 
are summarised in Table 1. We then fitted the intrinsic 
fluxes of the whole PWN with a PL F v oc v~ a . Like for 
the pulsar, we computed the fitting with a least-squares 
fitting method that takes into account the uncertainties of 
the measured fluxes as their weights. This yields a PWN 
spectral index of a™j fl = 0.56 ± 0.03 for the data in the 
nIR/optical range. 

We also qualitatively searched for evidence of possible 
variations in the nIR spectrum of the PWN as a function 
of position on angular scales smaller than ~ 1". We note 
that a systematic work done in the optical using HST data 
shows that the PL spectrum of the PWN does not change 
significantly when different positions are sampled both 
around the bar and around the putative jet and counter- 
jets (see, e.g. Serafimovich et al. 2004; Lundqvist et al. 

2011) . Fig.[5]shows the composite JHKs-band image of the 



Table 1. Extinction-corrected fluxes for the PSR 
B0540-69 PWN. Fluxes have been measured for the VLT 
J, H and Ks bands (first three rows) and for the HST* 
data presented in Mignani et al. (2010). 



Telescope/ 


Band 


Flux 


Instrument 




[1CT 27 erg s _1 cm" 2 Hz" 1 ] 


VLT/ NACO 


J 


5.41 ±0.22 




H 


6.75 ±0.27 




K s 


8.01 ±0.35 


HST/ WFPC2 


F336W 


2.11 ±0.22 




F450W 


3.61 ± 0.22 




F555W 


3.46 ±0.17 




F675W 


4.72 ±0.18 




F814W 


4.81 ±0.14 




F791W 


4.26 ±0.13 



*A11 the HST data were acquired in 2007, except for the data 
in the F791W filter that were taken in 1999. 

PSR B0540— 69 PWN, where the emission contributions 
in different bands are colour-coded. As seen, we cannot 
find evidence of regions more prominent in a band than in 
another. Thus, we conclude that there are no significant 
small-scale variations in the nIR spectrum of the PWN, 
like in the optical. We note, however, that the intrinsic 
variability of the PWN on angular scales of <^ 1" (e.g., De 
Luca et al. 2007) can affect the local multi-band spectral 
analysis. For this reason, the comparison of multi-band 
observations taken years apart must be taken with due 
caution. Indeed, while our NACO data have been taken 
within few weeks from each other (Sectn. 2.1), the avail- 
able HST data are about 3.5 to 15 years older (Mignani 
et al. 2010) and, in principle, are not directly compara- 
ble with the former. In the X-rays, a spectral analysis of 
different regions of the PWN has been performed using 
Chandra data (Petre et al. 2007). A quite complete dis- 
cussion of the comparison between the spatially-resolved 
optical and X-ray spectra of the PWN is given both in 
Serafimovich et al. (2004) and Lundqvist et al. (2011), 
and we will not repeat it here. 

4. Discussion 

4.1. The pulsar infrared luminosity 

We computed the PSR B0540— 69 luminosity by us- 
ing as a reference the most recent value of the LMC 
distance (48.97 ± 0.9 kpc), obtained from re-calibrating 
the Cepheids period- luminosity relation (Storm et al. 
2011). At the assumed distance, the measured extinction- 
corrected K-band flux of PSR B0540— 69 corresponds to 
a luminosity L K ~ 2.8 x 10 33 erg s" 1 . We plotted in Fig.[B] 
the K-band luminosity of all rotation-powered pulsars as 
a functional of the rotational energy loss rate E. We have 
corrected the observed K-band fluxes using the interstel- 
lar reddening Ay measured along the line of sight (see 
Mignani et al. 2007a and references therein) and the inter- 
stellar extinction coefficients of Fitzpatrick (1999). We as- 
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Fig. 5. Co-added JHK s -band image of the PSR 
B0540— 69 nebula. Images in different bands have been 
colour-coded (J=blue, H=green, Ks=red). 

sumed, the best available measurements of the pulsar dis- 
tances, either from optical/radio parallax or from the DM. 
We also plotted the measured upper limits for other pul- 
sars observed in the nIR, i.e. PSR J1119-6127 (Mignani 
et al. 2007a), PSR J1811-1925 and PSR J1930+1852 
(Kaplan & Moon 2006). As seen, the K-band luminosity 
correlates with the E, with Lk oc e 1 - 72 ^- 03 . The avail- 
able luminosity upper limits are also consistent with such a 
trend. The slope of the Lk~E relation is somewhat steeper 
than previously measured by Mignani et al. (2007a), i.e. 
L K oc E 1 - 3±0M i when psR B0540-69 was yet to be de- 
tected in the nIR. The former is equal to the correspond- 
ing slope measured in the optical band for the same set of 
pulsars, where cc _gi-70±o.03_ Notwithstanding all the 
associated uncertainties and the very limited statistics, 
the observed correlation between the nIR luminosity and 
E clearly clearly suggests that the pulsar nIR emission is 
rotation-powered and rules that it is entirely produced by 
a debris disk. This emission may then be of either purely 
magnctosphcric origin or partly coming from a pulsar's 
wind termination shock region. 

As seen from Fig. H2 PSR B0540— 69 is the most lumi- 
nous pulsar in the nIR. For its rotational energy loss rate 
(E ~ 1.5 x 10 38 erg s _1 ), the derived nIR luminosity cor- 
responds to an emission efficiency r\ n iR = Lk/E k 10~ 5 . 
This is higher than the emission efficiencies computed 
for all the other pulsars detected in the nIR, for which 
77 ~ 10~ 9 -10~ 6 (sec also Mignani et al. 2007a). In partic- 
ular, PSR B0540— 69 turns out to be about an order of 
magnitude more efficient that its "twin", the Crab pulsar, 
whose location in Fig.[B]is slightly below the luminosity--© 



Fig. 6. Measured K-band luminosities of identified 
rotation-powered pulsars (filled circles) and available up- 
per limits plotted as a function of the rotational energy 
loss rate E (updated from Mignani et al. 2007b). Different 
pulsars are labelled accordingly (see legenda) . The plotted 
errors on the K-band luminosity account for the uncertain- 
ties on the flux measurements, the extinction correction, 
and the nominal pulsar distance. For PSR B0540— 69, the 
error is comparable with the size of the symbol. The star 
indicates the candidate nIR counterpart to the X-ray pul- 
sar IGR J18490-0000 (Ratti et al. 2010; Curran et al. 
2011). 



trend. Although the measured flux of PSR B0540— 69 is 
sensitive to the choice of the background area (see Sectn. 
3.1), it is clear that the difference in the inferred nIR 
emission efficiency is too large to be ascribed to an over- 
estimate of the pulsar flux due to an inaccurate subtrac- 
tion of the nebula background. We note that a larger emis- 
sion efficiency of PSR B0540— 69 with respect to the Crab 
is also observed in the optical (e.g., Zharikov et al. 2006). 
Thus, it has to be related to a difference in the intrinsic 
pulsar emission properties, both in the nIR and in the op- 
tical. Non-thermal activity of young pulsars like the Crab 
and PSR B0540— 69 is due to spatially extended outer 
gaps in their magnetospheres. Unfortunately, there ex- 
ists no well-established self-consistent model of outer gaps 
which may be applicable even to pulsars of similar P and 
P. It is then difficult to explain, for instance, the difference 
between PSR B0540— 69 and the Crab pulsar in 7-rays. 
The same remains true for the optical to nIR emission. As 
seen from Fig.|51 the Vela pulsar, is clearly under-luminous 
with respect to the luminosity-© trend, which implies an 
intrinsically lower emission efficiency r\ n iR ~ 10~ 9 . A sim- 
ilarly low emission efficiency is also observed in the opti- 
cal (Zharikov et al. 2006). Like for the Crab, this must 
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be related to a difference in the intrinsic pulsar emission 
properties in these spectral ranges. 

We also compared the inferred nIR luminosity of 
the candidate counterpart to the X-ray pulsar IGR 
J18490-0000 (Ratti et al. 2010; Curran et al. 2011) with 
those of the pulsars identified in the nIR. As a reference to 
compute the IGR J18490— 0000 nIR luminosity, we used 
its observed magnitude Ks = 17.2 ± 0.4 and interstellar 
reddening E(B — V) — 6.62 (Curran et al. 2011) and we as- 
sumed a distance of 7 kpc, as in Gotthelf et al. (2011). The 
location of the IGR J18490-0000 candidate nIR coun- 
terpart in Fig. [S] is well above the luminosity- E trend, 
while its inferred nIR luminosity is comparable to that of 
PSR B0540-69 but for a factor of 10 lower E (Gotthelf et 
al. 2011). This would imply a nIR emission efficiency for 
IGR J18490-0000 which is a factor of 10 larger than PSR 
B0540— 69. Such a value is surprisingly high, even account- 
ing for the large uncertainties on the distance and inter- 
stellar extinction towards IGR J18490-0000 (Gotthelf et 
al. 2011). The inferred efficiency would be even higher if 
one assumes the observed magnitude (Ks = 16.4 ± 0.1) of 
the candidate counterpart originally reported in Ratti et 
al. (2010). This might argue against the proposed identi- 
fication, suggesting that the candidate nIR counterpart is 
an unrelated field star, which is plausible according to the 
relatively high (~ 6%) chance coincidence probability in 
the field (Curran et al. 2011). New nIR observations are 
required to obtain a more precise estimate of the distance 
and the extinction towards IGR J18490— 0000 and deter- 
mine the spectrum of its candidate counterpart. 

4.2. The pulsar multi-wavelength spectrum 

The multi-wavelength spectral energy distribution (SED) 
of PSR B0540-69 is shown in Figffl where the optical- 
to-nlR PL is compared with that best-fitting the Chandra 
spectrum (Kaaret et al. 2001) and its extrapolation in 
the optical domain. Due to the agreement between the 
optical and nIR spectra, the overall SED picture does 
not change with respect to that described in Mignani 
et al. (2010), with the spectral index of the optical-to- 
nlR PL (ao,niR = 0.70 ± 0.04) being quite similar to the 
X-ray one obtained from phase-integrated spectroscopy 
(ax = 0.92 ± 0.11), and the optical-to-nIR spectrum un- 
dershooting by a factor of 100 the low-energy extrapola- 
tion of the 0.1-10 keV X-ray spectrum. The peculiarity of 
the PSR B0540-69 SED, with the presence of two spec- 
tral breaks, has been already noted by Serafimovich et 
al. (2004) and Mignani et al. (2010), and we refer to the 
discussion therein. Here, we compare the optical-to-nIR 
spectrum of PSR B0540— 69 with those of other pulsars 
detected in both wavelength ranges. 

Although the correlation with the E (Sectn. 4.1) proves 
that the pulsar nIR emission is rotation-powered, all pul- 

3 A much lower probability (~ 1.6x 10 -5 ), however, is quoted 
by Ratti et al. (2010), although they do not explain how they 
derived this value. 




10" 10 15 10 16 10 17 10 1B 10 18 



v[Hz] 

Fig. 7. Optical-to-nIR spectral energy distribution of PSR 
B0540— 69 (points) compared with the PL model (Kaaret 
et al. 2001) best-fitting the Chandra X-ray spectrum (solid 
line) and its extrapolation in the optical domain. Dashed 
lines correspond to a la uncertainty on the model param- 
eters. 



sars feature markedly different behaviours in the nIR/mIR 
with respect to the optical/nUV, which makes it impos- 
sible to fit the overall spectrum with a single-component 
model. For instance, the Crab pulsar shows a clear spec- 
tral turnover red-ward of the nUV, where the PL spec- 
tral index a n uv gradually turns from positive into neg- 
ative, smoothly joining the optical/nIR spectrum (see, 
e.g. Fig. 4 of Mignani et al. 2010), but features a new 
change in the mIR, where it becomes approximately zero 
(Tcmim et al. 2009). For the Vela pulsar, instead, a sin- 
gle PL fits the flat nUV-to-nIR spectrum (Kargaltscv & 
Pavlov 2007) but the PL spectral index turns into posi- 
tive at wavelengths longer than 2.2 /zm, producing a spec- 
tacular rise of the PL spectrum in the mIR (Danilenko 
et al. 2011). For PSR B1509-58, the lack of multi-band 
detections at optical wavelengths, with only an R-band 
flux measurement obtained from polarimetry observations 
and, probably, affected by the presence of a brighter, 
nearby, star (Wagner & Seyfert 2000), does not allow to 
draw any definite conclusion. Nonetheless, the existing 
flux measurements are consistent with a single PL with 
positive spectral index fitting the optical/nIR spectrum. 
On the other hand, for PSR B0656+14 the composite op- 
tical/nUV spectrum sharply turns into a PL with positive 
spectral index in the nIR (Koptsevich et al. 2001), while 
for Geminga it shows a more gradual but similar turnover 
(Koptsevich et al. 2001) which seems to smoothly extend 
to the mIR (Danilenko et al. 2011). However, the cases of 
PSR B0656+14 and Geminga are different from those of 
the younger pulsar, since the observed spectral turnover 
between the optical/nUV and the nIR purely marks the 
transition between a Rayleigh- Jeans and a synchrotron- 
dominated continuum, and not between different slopes 
of the synchrotron PL spectra. 
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Thus, the spectra of the Crab and Vela pulsars sug- 
gest that PL spectral breaks between the optical/nIR and 
the mIR can be a characteristic of the pulsar magneto- 
spheric emission. This might be true, in principle, also 
for PSR B0540-69 and PSR B1509-58 for which, how- 
ever, no mIR observations are available yet. Similar breaks 
are also observed between the optical/UV and the X-ray 
bands (Mignani et al. 2007b) and are likely related to the 
different energy and density distribution of particles in 
different regions of the pulsar magnetosphere, where the 
optical and X-ray radiation is produced. The presence of 
spectral breaks also between the optical/nIR and the mIR 
suggests rather different emission regions for the observed 
radiation, while, on the other hand, the continuity be- 
tween the optical and nIR suggests that it is produced in 
very close regions in the neutron star magnetosphere. The 
comparison of the pulsar light curves in the optical and 
in the nIR/mIR would be crucial to verify this interpreta- 
tion, since differences in the light curve morphology and 
in the peak alignment would automatically reflect differ- 
ent locations of the emission regions. However, only very 
few pulsars have been observed to pulsate in the opti- 
cal (see Mignani 2010, 2011b for recent reviews), includ- 
ing the Crab and Vela pulsars and PSR B0540-69. Of 
them, only the Crab (Eikenberry et al. 1997) has been 
observed to pulsate in the nIR, although not yet with 
newest-generation instruments. For the Crab, the pulsar 
light curve evolves from the optical to the nIR, with evi- 
dence of an increase of the peak-to-peak phase separation 
and width as a function of wavelength, suggesting a differ- 
ent geometry of the emission regions. No pulsar has been 
searched yet for pulsations in the mIR. 

The level and spectral shape of intrinsic magneto- 
spheric emission in outer gap models depends strongly on 
many factors which may vary from pulsar to pulsar. In par- 
ticular, the details of the gap structure and its boundary 
conditions may be different, leading in turn to quite differ- 
ent properties of the non-thermal pulsar radiation. This is 
equally true in the high-energy domain (from X-rays to 7- 
rays) and in the low-energy domain (from UV to optical, 
to nIR and to longer wavelengths) . The most sophisticated 
3D model to date has been worked out for the Crab pul- 
sar only (Hirotani 2008). However, even this model has too 
many free parameters and, therefore, is subject to continu- 
ous modifications driven by new discoveries (Section 8.2 in 
Aleksic et al. 2011). Any high-quality photometric and/or 
spectral data in the low-energy range may shed light on 
the location and spatial extent of the regions of origin of 
emission which in outer-gap models comes from high-order 
generation cascades inside and outside the gaps. 

4.3. The PWN multi-wavelength spectrum 

Our NACO detection of the PWN around PSR B0540-69 
makes it one of the very few PWN systems detected both 
in the optical and either in the nIR or mIR. The others 
are the Crab (Sandberg & Sollerman 2009; Temim et al. 



2009), PSR J0205+6449 (Shibanov et al. 2008; Slane et 
al. 2008), and PSR J1124-5916 (Zharikov et al. 2008; 
Zyuzin et al. 2009). No other PWNe are detected so far 
in either of these spectral ranges, with the only exception 
of the PSR J1833-1034 PWN, which has been detected 
both in the nIR and mIR with the VLT and .Spitzer, re- 
spectively (Zajczyk et al. 2011). In all cases, the PWNe 
are also detected in the X-rays, where they have simi- 
lar morphologies as in the optical/IR, although, in some 
cases, with a somewhat larger spatial extent (e.g., the 
PSR J0205+6449 PWN). The relation between the op- 
tical/IR and X-ray spectra of these PWN, however, is 
not unique. For instance, for both the Crab and the PSR 
J0205+6449 PWNe, the nebula spectrum is compatible 
with a single PL, extending from the X-rays to the op- 
tical/IR. On the other hand, for the PSR J1124-5916 
PWN, the optical/IR spectrum clearly undershoots the 
extrapolation of the X-ray one, implying a break in the 
multi- wavelength PL spectrum between the optical/IR 
and the X-rays bands. 

Fig. |U shows the SED of the PSR B0540-69 PWN. 
Interestingly, the optical/nIR PL spectrum (a^nm = 
0.56±0.03) is flatter than the pulsar (a ,„m = 0.70±0.04), 
suggesting that the emission becomes harder at large dis- 
tances from the pulsar. We note that the PWN spectral 
index in the optical/nIR is flatter than measured in the 
optical by Serafimovich et al. (2004) based on a subset of 
our ffiTdata set. This is probably due to the 336W flux, 
which is lower than the best-fit optical/nIR PL and, then, 
produces a steepening of the PL in the optical. The best- fit 
Chandra PL spectrum of the PWN in 0.6-10 keV energy 
range Kaaret et al. (2001) is also shown in FigJSl together 
with the mIR fluxes from the Spitzer IRAC and MIPS ob- 
servations (Williams et al. 2008) and the radio fluxes from 
Manchester et al. (1993). The extrapolation of the PWN 
X-ray spectrum, with spectral index o^ WN = 1.04, to the 
nIR-optical domain clearly overestimates the PWN optical 
and nIR fluxes. This suggests that a spectral break must 
exist somewhere between the optical and X-ray bands. 
This behaviour is at contrast with what is observed in this 
wavelength range for PSR B0540— 69, where two spectral 
breaks are needed to match the pulsar PL spectrum be- 
tween the optical and X-rays (see Sect. 4.2). When extrap- 
olating the PWN nIR spectrum (ag^N = 0.56 ± 0.03) to 
the mIR, it is obvious that the fitted spectrum underes- 
timates the Spitzer fluxes. However, word of caution is 
necessary before drawing the conclusion that the spec- 
trum hardens between nIR and mIR bands. As described 
in Williams et al. (2008) the IRAC and MIPS mIR fluxes 
were estimated using rather large circular aperture (radius 
of 6") compared to the size of the elliptical aperture used 
in our calculations (see Fig(3]). Moreover, our nIR and the 
optical fluxes were corrected for the contribution from the 
pulsar and the background/foreground stars overlapping 
the PWN, while this is clearly not the case for the Spitzer 
fluxes, where both the pulsar and the stars are not re- 
solved at the angular scale of the IRAC and MIPS data. 
For this reason, the true intrinsic nebular fluxes in the mIR 
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Fig. 8. SED of the PSR B0540-69 PWN with our NACO JHK S fluxes (crosses), the HST optical fluxes (rectangles) 
obtained from the data of Mignani et al. (2010), the Spitzer IRAC and MIPS fluxes (triangles) of Williams et al. 
(2008), and the radio fluxes (diamonds) from Manchester et al. (1993). The Chandra PL fit to the X-ray spectrum 
of the PWN (Kaaret et al. 2001) is shown as a black bow-tie. The PL fit to the IR and optical data is shown as a 
dark grey bow-tie (see zoom on the right panel). The dotted line shows the extrapolation of the optical/nIR spectrum 
towards the X-ray and the mIR. 



are, most likely, lower that those presented in Williams et 
al. (2008). We note that the PWN spectral index has to 
change once more somewhere between the nIR/mIR and 
the radio bands from a^ N = 0.56 to o£ WN = 0.25 ± 0.1 
(Manchester et al. 1993). A similar spectral turnover to- 
wards radio frequencies is also observed for the Crab (e.g., 
review by Hester 2008; Serafimovich et al. 2004) and for 
the PSR J0205+6449 (Slane et al. 2008) PWNe. 

5. Summary and conclusions 

Thanks to high-spatial resolution observations performed 
with NACO at the VLT, we obtained the first detection 
of PSR B0540-69 and its PWN in the nIR. The pulsar's 
broad-band photometry fluxes in the optical (Mignani et 
al. 2010) and nIR are fit by a single PL with spectral in- 
dex ao,niR = 0.70 ±0.04, extending from ~ 0.3 to 2.8/xm. 
This implies continuity between the optical and the nIR 
spectra, although we cannot rule out yet the presence of 
a spectral break towards the mIR, like in the case of the 
Crab and Vela pulsars. The Ks-band luminosity of PSR 
B0540-69 (L K ~ 2.8 x 10 33 erg s" 1 ) well correlates with 
the rotational energy loss E, like in other pulsars, confirm- 
ing the non-thermal nature of the nIR emission. In partic- 
ular, at the LMC distance, PSR B0540-69 is the pulsar 
with the largest nIR luminosity, about an order of mag- 
nitude larger than the Crab's. This implies a correspond- 
ingly higher emission efficiency, like in the optical (e.g., 



Zharikov et al. 2006). The PWN has both a nIR morphol- 
ogy and a spatial extent consistent with those observed in 
the optical with the HST. We could not detect the bright 
emission "knot" observed in the HST images SW of the 
pulsar (De Luca et al. 2007). This can be explained ci- 
ther by a sharp break in its spectrum or, more likely, by 
its intrinsic variability, as shown by the multi-epoch HST 
observations (De Luca et al. 2007; Mignani et al. 2010). 
The optical/nIR spectrum of the PWN is also modelled 
with a PL (ag^N = 0.56 ± 0.03) which undershoots the 
extrapolation of the X-ray spectrum. This trend, so far 
observed only for the PSR J1124-5916 PWN, implies a 
single break in the optical/IR-to-X-ray spectrum of the 
PSR B0540-69 PWN, at variance with the pulsar for 
which two breaks are required. The optical/nIR spectrum 
of the PWN is also slightly flatter than the pulsar's, sug- 
gesting a possible spectral hardening at larger distance 
from the latter. Like in the optical, we could not find any 
obvious evidence of nIR spectral variations in the PWN 
on angular scales smaller than ~ 1". Only future observa- 
tions with the JWST will make it possible to extend the 
study of PSR B0540-69 and its PWN towards the mIR 
and search for a break in the pulsar spectrum at wave- 
lengths up to 5/im, possibly establishing a similarity with 
the Crab and Vela pulsars. At the other end of the op- 
tical spectrum, multi-band observations in the nUV with 
the HST, never performed so far, are required to search for 
the expected turnover between the optical and the X-rays. 
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